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g i 1. INTRODUCTION 

Superconductivity and Fermi superfluidity have been subjects of increasing in- 
terest since their discovery It seems as if nature wisely entangled Bose-Einstein (BE) 
condensation with these collective phenomena which are direct macroscopic manifes- 
tations of quantum physics. Recently, the nature of Cooper pairs in the BCS-BEC 
crossover has regained attention due to the observation of a large fraction of pre- 
formed fermion pairs on the BCS side of the Feshbach resonance in ultracold atomic 
Fermi gases [1]. While several theoretical explanations were proposed [2], the in- 
terpretations are still controversial. The root of the controversy is understanding 
what represents a Cooper pair in a correlated Fermi system, a concept that has been 
extensively used in the literature to describe the properties and dynamics of the 
superconducting state but without a unique and precise definition. 

Although the question What is a Cooper pair? seems academic, it is indeed a 
question of upmost practical importance given the great advances we are witnessing 
in ultracold atom physics and quantum information. New experimental probes allow 
access to unusual regimes of strongly-coupled systems with unprecedented degree of 
control. Therefore, concepts need to be sharply defined to correctly interpret these 
new experimental results. The truth is that much of our current physical intuition is 
built upon mean-field solutions. Can one get new useful insight out of exactly-solvable 
models? This manuscript discusses the Cooper-pair concept at the very fundamental 
level and proposes a qualitatively different scenario of the BCS-BEC crossover, based 
on the exact solution to the BCS Hamiltonian [3]. Only in the extreme BEC limit 
does this new scenario and Leggett's naive ansatz [4] become identical. As we will 
see, the Cooper pair in a correlated Fermi system is merely a statistical concept 
whose manifestation and interpretation depend upon the particular measurement or 
physical operation performed in the system. 



While the superconducting and Fermi superfluid thermodynamic states represent 
a mixed-system of quasifree and pair-correlated fermions, the molecular BEC which 
arises in the dilute and strong coupling limit has all fermions bound into pairs forming 
a unique macroscopic quantum state. It is by now well accepted in which sense these 
states represent a Bose-Einstein (BE) condensation. What defines a BE condensa- 
tion in an interacting iV-particle system is spontaneous gauge symmetry breaking, 
or phase coherence, of its quantum state (a concept that strictly applies in the ther- 
modynamic limit (TL)). Yang [5] provided a detailed mathematical characterization 
of this phenomenon which manifests itself as off-diagonal long-range order (ODLRO) 
or, equivalently, by the existence of an eigenvalue of order N in a reduced density 
matrix p n , defined as p n (kl • • • \ij • ■ •) = Tr[a,k(ii • • • p ■ ■ ■ a-a\\ (i,j, k, /, • • •, labelling 
single-particle states annihilated by bosonic and/or fermionic operators a's) where 
n is the number of particles forming the smallest unit or complex that condenses. 
(n = 2 electrons for BCS superconductors). 



2. WHAT IS THE MANY-BODY WAVEFUNCTION? 

The above definition of a BEC does not imply the naive view of a many-body 
state of quantum objects with identical internal wavefunctions. Indeed, we will argue 
that a current-carrying ground state (GS) of a uniform superconducting or Fermi 
superfluid iV-particle system is of the form 

ty(xi,---,x N ) -A [4>i(xi,x 2 ) ■■■4>n/2(xn-i,xn)] , (1) 

with Xj = (tj, o~j), antisymmetrizer A, and the pair state 

• ( r i+ r j) 

where x is a spin function of a certain symmetry, q is the pair center-of-mass mo- 
mentum, and (p a {r) the internal wavefunction which may represent either a quasi- 
molecular resonant state or a scattering state (i.e., a mixed state), depending upon 
the strength of the interaction between particles. The antisymmetrizer A is applied 
to the particle index x^. For example, \l/(xi,X2, • • • , xn) = —^(x2, %i, • • • , %n)- 

For pedagogical reasons we will recall the one-Cooper-pair problem and then 
address the question of what happens when we add more pairs (q = in the GS). 
The Cooper-pair solution can be recovered from the wavefunction (1) by assuming 
that N — 2 fermions cj^ are free, filling the lowest momentum states k (of energy e^) 
up to the Fermi momentum kp, while only an additional fermion pair (with momenta 
k > &f) feels an attractive s-wave interaction in the spin singlet channel 
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The role of the Fermi sea, \F) = rifc<fc F c kT c -ki 1^)' ^ s ^° Pauli-block states below the 
Fermi energy ep. Assuming that the attractive pairing interaction G < is constant 
around the Fermi energy, the eigenvalue E satisfies 
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which always has a real negative solution, implying that the Cooper pair is bound 
regardless of the strength of the attractive interaction. The Fermi sea is then unstable 
against the formation of bound pairs of electrons. 

What happens when the pairing interaction also affects electrons in the Fermi 
sea? The answer to this question is the BCS theory whose canonical form considers 
a simple variational state of the form (1) with all identical internal wavefunctions (a 
pair mean-field) 
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where Pm is the projector onto the state with M fermion pairs, and i>k, u^ are the 
BCS amplitudes 
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with A the gap and (i the chemical potential. In the language of first quantization 
this mean-field state represents a Pfaffian state: Given the skew-symmetric matrix 
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(xn-i,xn)], which represents a linear 



Pf\& = *bcs(^i, • • • , %n) = A [0(xi, x 2 ) 
combination of (N — 1)!! terms. 

Since the structure of the BCS pair is averaged over the whole system, it does 
not characterize a Cooper pair in the superconducting/superfluid medium except in 
the extreme strong-coupling and dilute limits where all pairs are bounded and non- 
overlapping. (This approximate state is not the exact ground state of the reduced 
BCS Hamiltonian, Eq. (3).) Sometimes, the structure of the Cooper pair is associated 
with the pair-correlation function (BCS|cL.cl k , |BCS) = u^v^ describing the pair 
correlation among fermions of the same pair as well as the exchange between fermions 
of different pairs. 

What is the nature of a Cooper pair in a Fermi superfluid or superconducting 
medium? To address this question we will use the exact solution to the reduced BCS 
Hamiltonian 

Q 

H = y^ £ k nu + — Y^ cLc' uc_ t /ic w , (3) 
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proposed by Richardson 40 years ago [3,6]. H involves all terms with time-reversed 
pairs (k f, — k |) from a contact interaction. It is consistent with an effective single- 
channel description of the BCS-BEC crossover theory [4] in terms of a zero-range po- 
tential. Realistic finite-range interactions produce qualitatively similar results along 



the crossover [7]. The main differences are expected in the BEC region where the re- 
duced BCS Hamiltonian, Eq. (3), cannot account for the reduction in the scattering 
length of the composite pairs [8] and the corresponding depletion of the molecular 
condensate [9]. 

For simplicity we will consider N = N^ + Ni spin-1/2 (i.e., 2-flavor) fermions in 
a three-dimensional box of volume V with periodic boundary conditions, interacting 
through an attractive constant (s- wave- singlet-pairing) potential. (Other pairing 
symmetries can also be accommodated.) Exact N = 2M + v particle eigenstates of 
H can be written as 

M 

m = n s » > with s t = H^—w c ^ c] -^ (4) 

1=1 k k £ 

where \u) = \v±, v^ ■ • ■ , vl) is a state of v unpaired fermions [y = J^ k v^ with u^ = 
1, 0) defined by c_ k ,c^ \v) = 0, and n^ \v) = u^ \u). L is the total number of single 
particle states. The GS |^ ) is m the v = (iV|=iVjJ sector. 

Each eigenstate |\&) is completely defined by a set of M (in general, complex) 
spectral parameters (pair energies) Eg which are a solution of the Richardson's equa- 
tions 

i + -y^^ + — y — - — = o, (5) 

k rn(^t£) = l 

and the eigenvalues of H are S = J^ k £k^k + J2i=i Ee [6,10]. A crucial observation is 
that if a complex E(, satisfies (5), its complex-conjugate E\ is also a solution. Thus, 
|\1/) restores time- reversal invariance. The ansatz (4) is a natural generalization of the 
Cooper-pair problem without an inert Fermi sea, and with all pairs subjected to the 
pairing interaction. The pair structure in (4) is similar to the Cooper pair in (2), and 
the many-body state has the form (1) with all pairs different. In a sense, that will 
become clear later on, the many-body state (4) represents a statistical distribution 
of resonances (Schafroth's intuition of the superconducting state). 



3. THERMODYNAMIC LIMIT OF RICHARDSON'S EQUATIONS 

Since we are concerned with uniform bulk Fermi systems, we are interested in the 
TL (i.e., N, V — > oo with p = N/V = kp/dn 2 = constant) of the Richardson's exact 
solution Eq. (5). This limit was studied by Gaudin [11] when the energy spectrum 
O G [—a;, a;] is bounded, and v = 0. Eqs. (5) reduce to the well-known gap equation 

\tde , 9{£) +1=0, (6) 



2jn yj{e - M ) 2 + A 2 G 

where g(e) represents the density of states. In our case Q G [0, oo] is unbounded 
with g(s) = m 3 / 2 ^ 1 / 2 /(\/2n 2 h ) for e^ = % k 2 /2m. Due to the absence of an upper 
cutoff, Eq. (6) is singular and the TL becomes a subtle mathematical procedure 
whose solution will be presented here. 



There are two ways to regularize the problem: One can either introduce an energy 
cutoff in the integrals or one can cancel the singularities with physical quantities 
whose bare counterpart diverges in the same way [4]. For this problem, Eq. (3), the 
physical quantity is the scattering length a s given by 
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The non-singular gap equation (after integration [13]) is 
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where energies are now in units of ep = h k F /2m and lengths in units of £p = 1/kp. 
Pp{x) is the Legendre function of the first kind of degree (3. The equations for the 
conservation of the number of pairs M 



4 

3^ 



W + ^ 2 + A 2 ) s / 4 P 



and GS energy density [12], for arbitrary strength r], 
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do not need regularization, these are exact results. Indeed, for a given rj, the chemical 
potential \i and gap A need to be determined self-consistently from Eqs. (7) and (8) 
(Fig. 1). Then, the GS energy can be computed as a function of density using Eq. 
(9) . It shows no phase segregation. The exact binding (or condensation) energy per 
electron pair Eb is 
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One can analytically determine the asymptotic rj — ► ±oo behavior of \i and A 
by properly combining Eqs. (7) and (8), and making use of the defining Eqs. for the 
Legendre functions Pp{x). Indeed, one can convert Eqs. (7) and (8) into coupled 
differential equations, with one particular combination leading to 



d\n(ix 2 + A< 



7rA 



!+(¥)£ 



The asymptotics can be determined by using this equation and the results are de- 
picted in Table 1. 

The complete solution of Eqs. (5) in the TL [11] amounts to determining (for a 
given rj) the set of pair energies Eg which, for the GS, form a single open, continuous, 



Table 1 



Analytic expressions for selected values of r\ = l//cpa s ; x is the root of Pi, i.e. 
Pi(x) = 0, and E(y) is the complete elliptic integral of the second kind. Note that 
5n 2 £ (r] = 0)/k% = /j(r] = 0) w 0.590606. 
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and reflection-symmetric arc F with extreme points Ep = 2(p±iA). Here, we simply 
present the results without derivation. The equation for F is 
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where z = y (E — jj) + A 2 . 



4. BCS-TO-BEC CROSSOVER PHYSICS 



Figure 1 shows the BCS-BEC crossover diagram and several arcs corresponding 
to particular values of rj. Since a crossover diagram is not a phase diagram, i.e. 
there is no symmetry order parameter or non-analyticities sharply differentiating the 
regions, boundaries are in principle arbitrary. 

Here we have adopted the following criteria: The geometry of the arcs F serves 
us to establish a criterium for defining boundary regions in the crossover diagram. 
In the extreme weak-coupling limit, rj — ► — oo, the pair energies are twice the single 
particle energies (Eg — ► 2£k). As soon as the interaction is switched on, i.e. kpa s is 
an infinitesimal negative number, a fraction f of the pair energies close to 2/j, become 
complex, forming an arc F in the complex plane. The fraction 1 — f of fermion pairs 
below the crossing energy 2e c of F with the real axis have real pair energies. They 
are still in a sea of uncorrelated pairs with an effective Fermi energy e c . The dark 
grey region labelled BCS, which extends from rj = — oo to r] = 0.37, is characterized 
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Figure 1. Ground state crossover diagram displaying the different regimes as a 
function of r\. A few arcs T, from Eq. (10), corresponding to the values of rj = -1, 0, 
0.37, 0.55, 1 and 2, whose extremes are Ep = 2(/t ± iA), are displayed. The dashed 
line corresponds to the Feshbach resonance r\ = 0. Open circles represent the values 
of E for which Cooper-pair wavefunctions are plotted in Fig. 3. Insets show the 
behavior of the gap A and chemical potential \x. 



by a mixture of complex pair energies with a positive real part and real pair energies. 
rj = 0.37 is the value at which all pair energies are complex, i.e. f = 1, and the 
effective Fermi sea has disappeared (e c = 0). Within the BCS region we plotted the 
arcs for rj = —1 (solid line), with a fraction f = 0.35 and an effective Fermi energy 
e c = 0.747 , and rj = (dashed line) with f = 0.87 and e c = 0.25. The pseudogap 
region P, indicated in light grey, extends from rj = 0.37 to rj = 0.55 where \x = 0. 
Within this region the real part of the pair energies changes from positive to negative, 
and P describes a mixture of Cooper resonances and quasibound molecules. The BEC 
(white) region, rj > 0.55, is characterized by all pair energies having negative real 
parts, i.e. all pairs are quasibound molecules. Within the BEC region we plotted 
two arcs with rj = 1 and rj — 2. As r\ increases further, F tends to an almost vertical 
line with Re(E) ~ 2^, and -2A < Im(E) < 2A. 



5. WHAT IS THE STRUCTURE OF A COOPER PAIR? 



Having established a qualitative BCS-BEC crossover diagram, we turn now to 
the nature of the Cooper pairs. If the literature is not clear about the size £ of a 
Cooper pair (all lenghts are written in units of £f = 1//cf), it is partly because 
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Figure 2. Different definitions of Cooper-pair sizes £, and condensate fraction A and 
f as functions of the interaction strength r\. The horizontal solid line represents the 
interparticle distance r s = (97T/4) 1 / 3 (lenghts are written in units of ^p = l/kp). £e 
and £p increase for negative r\ values, and eventually diverge when r\ — > — oo as the 
inverse of the gap, like Pippard's coherence length £ — 2/(ttA). 



it is not clear what a Cooper pair is. Pippard, in his nonlocal electrodynamics 
of the superconducting state, introduced the characteristic length £o by using an 
uncertainty-principle argument involving only electrons within a shell of width 2A 
around the Fermi energy. The resulting quantity, named Pippard's coherence length, 
£o = 2/(7rA), is usually associated to £. On the other hand, an acceptable definition 
could be £ = \J {(p\r 2 \(p) , (<p\<p) = 1. From Eq. (4), the Cooper-pair wavefunction is 
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constant. Thus 
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for arbitrary interaction strength r\. In the weak-coupling BCS limit (A <C \i rj 1), 
when E = Ep, we get £e = tt^o/ v2. On the other hand, in the same limit, if one 
uses (^k = CpWkfk, one gets £p = £e/2, and if one uses (p^ cs = CBes^k/Mk, one gets 
£bcs = V21/2. 



Analytic expressions (using dimensional regularization) for £p and £bcs can also 
be obtained for arbitrary coupling strength. The results are 
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where 7 = — Pe(Pi?)/|Pi?|. Figure 2 displays the different sizes, all labelled by £, as a 
function of 77. Clearly, in (4) there is no unique pair size £e but a distribution of sizes. 
Although in the BCS region there is a distribution of sizes from the smallest pair with 
E = Ep in the extremes of the arc to an almost infinite size pairs when they are 
close to the real axis, already at resonance (77 = 0) most of the pairs have the same 
structure. This latter feature becomes more pronounced in the BEC region. Figure 
2 shows in solid line the smallest size in Y corresponding to Ep. On the other hand, 
there is a unique pair size for <pp and </?bcs; the dashed and dotted lines correspond 
to the £p and £bcs sizes respectively. Notice the relation between sizes and the 
interparticle distance r s = (97T/4) 1 ' 3 which is indicated as a full horizontal line. While 
£e and £p increase for negative 77 values, and eventually diverge when 77 — > —00, £bcs 
tends to a constant value, related to the interparticle distance, showing its essentially 
uncorrelated nature. The fact that £p < £_e/2 in the overlapping BCS region is 
a clear demonstration that £p measures the mean distance between time-reversed 
pairs irrespective of the Cooper pair they belong to. For non-overlapping pairs (BEC 
region) both sizes converge to the same values. 

The various definitions of Cooper-pair wavefunctions are depicted in Fig. 3, 
which compares </?BCs( r ) 5 Vp(?*), and ^Pe{ t ) for interaction strengths which corre- 
spond to the BCS, Feshbach resonance, and BEC regions of Fig. 1. Notice that 
while </?bcs? and <pp are highly oscillatory in the weak-coupling region, this is not the 
case with ipp which always has an exponential form. Clearly, a single and unified 
picture emerges for the crossover when using a many-body state such as (1): <pe 
evolves smoothly through the crossover as it should. It is important to mention that 
the three wavefunctions are exactly the same in the strong coupling limit 77 — > +00. 
In this limit the naive ansatz of Leggett [4] and the GS (1) coincide, becoming the 
simple Pfaffian state ^bcsO^I; ' ' ' ? x n) = A [4>{ x ii Xc i) ' ' '4>{ x n-i, x n)]- 



6. CONDENSATE FRACTION OF A FERMI SUPERFLUID 

The condensate fraction is a thermodynamic property charaterizing the super- 
fuid/superconducting state. In an ideal Bose gas it is the fraction of bosons occupying 
the lowest-energy quantum state. In a nonideal (Bose or Fermi) gas, one cannot use 
this definition. Following Yang [5], ODLRO in ^(ri j r 2 | | r 3 | r 4 |) — > (/>*(rit, r 2 | 
)</>(r 3 f , r 4 |) may be used to define 
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Figure 3. Cooper-pair wavefunctions for different 77. The upper and middle (Fes- 
hbach resonance) panels correspond to the BCS region, while the bottom one is in 
the BEC region. Except for the Cooper ifE case, the other two wavefunctions always 
vanish at r = 0. It is only in the limit 77 — > +00 that the three states exactly coincide. 



as a measurement of the condensate fraction in a two-flavor Fermi system. Figure 
2 shows A and the fraction f of (Cooper) pairs in the arc, that is, the fraction of 
complex spectral parameters. Although A can qualitatively describe the fraction of 
correlated pairs, the fraction f gives the value 1 at the BCS-pseudogap boundary 
(77 = 0.37) while A = 1 for 77 — > 00. We note that at resonance (77 = 0), f = 87% 
of the fermions form Cooper pairs (A ~ 0.7). These results provide a qualitative 
explanation of the experiments in [1] without resorting to a projection method [2]. 
Close to resonance on the BCS side, a fraction f ~ 80% of the atoms form Cooper 
pairs with sizes comparable to r s . Those atom pairs are efficiently transformed into 
quasimolecules during a rapid magnetic field ramping across the resonance giving rise 
to the molecular condensate fractions observed experimentally. 

Finally, Table 1 summarizes various analytical results for the chemical potential, 
gap, binding (or condensation) energy, and condensate fraction, in four important 
limits of the single-channel model: 77 — ► —00 (BCS), 77 = (Resonance), 77 = 
(BCS-to-BEC crossover), and 77 — > 00 (Extreme BEC limit). 



7. CONCLUSIONS 

In summary, we studied the BCS-BEC crossover problem, as well as the nature 
of Cooper pairs in a correlated Fermi system, from the exact GS \^o) of the reduced 
BCS Hamiltonian. We have analytically determined its exact TL for the quadratic 
single-particle dispersion, and calculated several quantities of physical interest, such 
as the binding (or condensation) energy and the condensate fraction. The funda- 
mental difference between the BCS state and \^o) is the fact that the former, being 
a mean-field, pairs off all electrons in a unique paired wavefunction while the ex- 
act state contains a statistical distribution of quasimolecular resonant and scattering 
states (i.e., a mixed state) that depends upon the strength of the interaction between 
particles. It is only in the extreme BEC regime that the two pictures coincide asymp- 
totically. Then, what a Cooper pair is depends upon the particular measurement or 
operation performed in the physical system. The validity of the present description is 
not restricted to integrable pairing Hamiltonians, but rather the ansatz \^o), which 
is a natural extension of the Cooper problem, could be used as a variational state for 
more realistic interactions. The Cooper-pair wavefunction f>E has a fixed functional 
s-wave form that accommodates pair-correlated fermions close to the Fermi energy 
in the weak coupling limit as well as quasibound molecules in the BEC limit. It 
is free from the characteristic oscillations displayed by v?bcs an d <£>p related to ex- 
change among pairs. The BCS region in the crossover diagram represents a mixture 
of Cooper pairs and quasifree fermions. These properties resemble the Schafroth de- 
scription of the superconducting state [16]. It may be argued that the single- channel 
model is insufficient to describe the system close to resonance where the degrees of 
freedom associated to the molecular closed channel should be explicitly incorporated. 
A derivation analogous to the one presented here can be pursued by using a recently 
proposed atom-molecule integrable model [14] which captures the essential features 
of the two-channel model [15]. The structure of the Cooper pair in this new model 
is the same as in (4). 
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APPENDIX: USEFUL INTEGRALS 

We simply quote the results of some integrals used in this article. Some of 
them were obtained using the technique of dimensional regularization (z = a + ib, 

z* = a — ib, 7 = —a/\z\, excluding Re(z) = a > and Im(z) = 6 = 0) 
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